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Abstract 

Potential energy surfaces calculated with self-consistent mean-field methods are a very powerful 
tool, since their solutions are, in theory, global minima of the non-constrained subspace. How- 
ever, this minimization leads to an incertitude concerning the saddle points, that can sometimes 
be no more saddle points in bigger constrained subspaces (fake saddle points), or can be missing 
on a trajectory (missing saddle points). These phenomena are the consequences of discontinuities 
of the self-consistent potential energy surfaces (SPES). These discontinuities may have important 
consequences, since they can for example hide the real height of an energy barrier, and avoid any 
use of a SPES for further dynamical calculations, barrier penetrability estimations, or trajectory 
predictions. Discontinuities are not related to the quality of the production of a SPES, since even 
a perfectly converged SPES with an ideally fine mesh can be discontinuous. In this paper we 
explain what are the discontinuities, their consequences, and their origins. We then propose a 
numerical method to detect and identify discontinuities on a given SPES, and finally we discuss 
what are the best ways to transform a discontinuous SPES into a continuous one. 
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Introduction 

Potential energy surfaces (PES) are a widely used tool to describe physical and chemical sys- 
tems, among others. For example, they are used in DNA chemistry yj, materials chemistry [^, 
chemical physics |3], astrophysics fl, etc. . . Numerous methods exist to extract from a PES local 
minima, saddle points, and least energy paths between local minima (EJ^. To produce a PES, a 
method giving the energy of a system as a function of a given number of constrained variables 
is needed. In this paper, we will separate these methods into two main classes: self-consistent 
and non-self-consistent methods. Surfaces associated with these methods will be called Self- 
consistent Potential Energy Surfaces (SPES) and Non Self-consistent Energy Surfaces (NSPES), 
respectively. In nuclear physics, SPES can be obtained for example by constrained Hartree- 
Fock methods and extensions [:3,SS[l2li by constrained Relativistic Mean Field method I fill , 
etc NSPES can be produced by several methods, ranging from the historical liquid drop 



model lll2Lll3[l to the well-known macroscopic-microscopic model, using parametrization of the 



nuclear mean-field deformation ifTil [isll . 
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The main difference between these two classes of methods is the way they deal with the 
non-constrained degrees of freedom of the system. Non self-consistent methods neglect their in- 
fluence (all non-constrained degrees of freedom take a fixed value), while self-consistent methods 
perform an automatic minimization of the energy of the system in the non-constrained subspace. 
The presence or the absence of this automatic minimization leads to one specific problem for 
each class of methods: 

• NSPES problem: every point of a given NSPES may not be a minimum of the same 
NSPES with an additional dimension. 

• SPES problem: a saddle point on a trajectory from a SPES is not automatically visible on 
the same trajectory from the same SPES restricted to a smaller dimension (missing saddle 
point), and a saddle point from a SPES is not automatically a saddle point on the same 
SPES with an additional dimension (fake saddle point). 

These two problems have important consequences. For the NSPES problem, one would 
expect that enriching a NSPES with an additional constraint while conserving a given set of 
symmetries should improve the quality of the description, instead of invaliding it, but this is not 
always the case. For example, quantities calculated on a A^-dimensional NSPES like minima, 
saddle points, or least energy paths between local minima may differ dramatically from quan- 
tities calculated on a (A^ -h l)-dimensional NSPES, with the first constraints and conserved 
symmetries being the same for both NSPES. Let us imagine an hypothetical non-constrained 
HFB minimum with (iAIQtoIiA) - 37 b^'^. Finding this minimum with non-self-consistent 
methods will require to explore the deformation degree of freedom corresponding to Qjq, with- 
out knowing beforehand that this degree of freedom has a role to play. If this degree of freedom 
is not taken into account, the resulting minimum can be anything but the searched minimum, and 
may differ dramatically from it. 

Concerning the SPES problem, the obtained saddle points can he fake or missing. This last 
problem is well known (it corresponds to a tipping point in Catastrophe Theory 1 16]), and is used 
sometimes to contest -with good reason- the quality of the results obtained by self-consistent or 



minimization methods BlTll . For example, in 111811 . the authors present a schematic case leading to 



the calculation of a fake saddle point by a minimization method, and conclude rather quickly that 
the doubt concerning the reality of saddle points found with self-consistent methods invalidates 
any production of a SPES by a minimization method: "[...] in Hartree-Fock-Bogoliubov (HFB) 
calculations with multiple constraints, saddle-point shapes and energies frequently cannot be 
determined accurately and are subject to errors of fairly random magnitude". 

To our knowledge, there is no numerical or analytical method to fix the NSPES problem. 
Any result extracted from a NSPES may be invalidated by a NSPES of higher dimension. By 
comparison, there exists a simple numerical method to fix the SPES problem, that we will present 
in this paper 

A saddle point that is no more a saddle point in a bigger constrained subspace and a saddle 
point which vanishes while reducing the constrained subspace are two of the possible conse- 
quences of a unique phenomenon that we call a discontinuity of the SPES. A continuous SPES 
has no fake or missing saddle point. In this paper we first show such discontinuities, explain some 
of their consequences, and explain their origins. We then propose a method to quantitatively es- 
timate if a SPES is continuous or not, and to identify the discontinuities if present. Finally we 
discuss what are the best ways to get a continuous SPES from a discontinuous one. 
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1. Discontinuities 



1.1. Model used 

Every calculation presented in this paper has been obtained by a two-center basis Hartree- 
Fock-Bogoliubov model |19] with Gogny DIS nucleon-nucleon effective interaction 211. 



based on the minimization principle of the energy functional, namely 

5{<P\H - AnN -AzZ-J^ AiQ,o\ilr) = 0, (1) 

where H is the nuclear microscopic Hamiltonian, Qio is a mass multipole operator, and A^, 
Az and A/ are the Lagrange multipliers corresponding to constraints on the neutron number A^, 
proton number Z, and mass multipole moment qio, respectively. The corresponding constrained 
equations are 

((A|Z|<A> = Z, (2) 
(ifrim) = A^, (3) 
(-Aiaol^) = qio, (4) 

with Qio defined as 

Qio = (I + 61,2)^1^^ J]r'.Y,o(0i,<f>i)- (5) 



!=1 

The HFB energy of the constrained system \iff(N,z.{qio})) is 

^fN,l[q,o\) ^ {^(N,Z,{qM\H\ilf(N,Z,{qia)))- (6) 

The HFB equations (Eqs. ([T]i to (|4|i) are solved iteratively by expanding the single-particle wave- 
functions onto an axial two-center harmonic oscillator basis. The basis parameters have been 
automatically optimized for each set of constraints. The Lagrange multipliers are adjusted at 
each iteration in order that the constraint conditions (Eqs. (|2]i to (HJi) are fulfilled, and the gener- 
alized density matrix is corrected accordingly. This technique allows us to solve the constrained 
HFB equations in the restricted variational space orthogonal to the PES for which the constraint 
conditions are fulfilled. In this way, the problems of non-uniform mesh and inaccessible re- 



gions mentioned in e.g. 112211 are completely avoided. For more information on the adjustment of 



the Lagrange multipliers in the case of multiple linear constraints, see for example Appendix A 



of II23I1. 



1.2. Example of a discontinuity 

During the production of a self-consistent potential energy surface, a strange behavior can 
sometimes be observed. For a given set of constraints, starting the iterative process from two 
different points can lead to two very diff'erent solutions, both converged. This is clearly not what 
is expected from a self-consistent method, which is supposed to always converge to the only 
solution minimizing the total binding energy under the action of the constraints, whatever the 
starting point. Such a strange behavior is shown on Fig. [1] where two one-dimensional SPES 
have been produced, each point being the starting point for its left (squares) or right (triangles) 
neighbor. One can see that for a rather large range of constraint values, two different solutions 
can exist. The "bow tie" shape formed by these two curves is characteristic of a discontinuity in 
the SPES. 
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1.3. Explanation 

For an hypothetical system, let us consider the two-dimensional SPES plotted as contour 
lines on the upper panels (a), (b) and (c) of Fig. |2l showing two valleys A and B in the (Qio, 
Qxx) subspace. Let us now try to produce the one-dimensional SPES with increasing values of 
^20, i e- each point being the starting point for the calculation of its right neighbor The resulting 
total binding energy is plotted on the panel (d), and the mean value {Qxx) is plotted as dots on 
the panel (a). One can see that the converged solutions follow valley A until the system pass the 
barrier and "falls" into valley B. If we do the same calculations with decreasing values of q2Q, 
we obtain the curves on panels (b) and (e). This time, the system stays longer in valley B before 
falling into valley A. If we superpose the two curves on the panels (d) and (e), we see the same 
"bow tie" shape as in Fig. [1] 

The typical hysteresis figure shown on the lower panel of Fig. [Uor in panels (a) and (b) of 
Fig.|2]is due to the fact that the system can stay in its valley of origin even when the solution is 
no longer the global minimum. This phenomenon is mainly related to the way the minimization 
is performed. In the ideal case of a hypothetical perfect minimization, the global minimum is 
always found, leading to curves in the panels (c) and (f). 

To solve the hysteresis problem in the case of a behavior like the one shown in Fig. [1] 
one has to discard the non-global minima, keeping the states with the lowest binding energies. 
This solution leads us to the previous ideal case: every constrained HFB solution is an absolute 
minimum, there is no more hysteresis, thus no more influence of the starting points. In the 
following, we will assume that the minimization is perfect, or that hysteresis problems have been 
tracked and solved. However, even with a perfect minimization, the passage from one valley to 
the other (from point a to point b) is not continuous, since the mean value {Qxx) is quite different 
for these two points. This discontinuity is not visible in the total binding energy curve, since 
point a and point b can be found to have very close (if not the same) total binding energies. A 
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Figure 2: (Color online) Illustration of the solutions found when the propagation is from the left neighbor (panels (a) and 
(d)), from the right neighbor (panels (b) and (e)), and in the ideal case of a perfect minimization (panels (c) and (f)). The 
upper panels show the same schematic two-dimensional SPES as contour lines, with two valleys labeled A and B. The 
lower panels show the total binding energy as a function of the Q20 values. In panel (f), the point with the maximum 
total binding energy can be the point a orb from panel (c). 

discontinuity in a SPES can have serious consequences. 
1.4. Consequences 

We define a path in a subspace as a 1 -dimensional set of contiguous points in this subspace. 
In the case of a non-regular mesh, the contiguousness condition becomes a closest neighbors 
condition (cf. 12. 3 1 for a definition of this condition). If a path in the constrained subspace crosses 
a discontinuity, there exists a subspace in which its set of points is not a path anymore, since the 
points before and after the discontinuity belong to different valleys and are no more contiguous. 
A description of the evolution of a system along a path that crosses a discontinuity misses the 
description of the passage between these valleys. This missing passage may correspond to a 
high-energy saddle point. 

One can generalize this statement to the case of A^-dimensional SPES: any description of the 
evolution of a system on a discontinuous SPES is an incomplete description, since when crossing 
a discontinuity, the real system does not "teleport" itself from one valley to the other 

An other consequence is that the energy barrier between valleys A and B is necessarily 
wrongly evaluated, since there can be a huge energy barrier or a low-energy path between both 
valleys. The barrier height that can be seen on the panel (f) of Fig.|2]may be totally different from 
the real barrier, and there is no way to estimate this difference by only considering the evolution 
of the total binding energy like on Fig.[T]or on the panel (f) of Fig.|2] In other words, the visible 
saddle point a or b on the panel (f) is a fake saddle point. If the evolution of the energy on the 
panel (f) were to be monotonic, a saddle point would not even be visible, and point a orb would 
correspond to a missing saddle point. 

2. Finding the discontinuities 

In this section we propose a method to estimate if a given SPES is continuous or not, and to 
localize and identify the discontinuities. 
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2.7. Definitions 

What is usually called a self-consistent potential energy surface (SPES) is a set of A' self- 
consistent solutions {\>l/i))i=\,N under the action of Nc constraint operators {Q.j}j=\,N, with con- 
straint values qi j. This set can be regularly or randomly distributed on the constrained subspace. 
We suppose that a boundary of this SPES is given, separating the deformation subspace in an 
interior 7 containing the SPES and an exterior. For example, a SPES with constraint values on 
220 regularly spaced from 100 b to 200 b with a step q^Q = 10 b has an interior defined as 
7= [100,200]. 

Let us introduce the idealized self-consistent potential energy surface (ISPES): for a given 
SPES, the unique corresponding ISPES is the infinite set of self-consistent solutions Wxjj'^)] under 
the action of the same Nc constraint operators [Q.j]j=\,N, with constraint values qxj taking all 
possible values in the same interior. The x quantity is a A^c-dimensional vector of coordinates. 
The ISPES corresponding to the SPES in the previous example can be obtained by taking the 
limit ~^ 0: 

= lim ({|^,)}). (7) 

We define the continuity of an ISPES: an ISPES is continuous if the mean values of every 
one-body multipole operator appUed to the system are continuous with respect to each constraint 
coordinate: 



We then define the continuity of a SPES: a SPES is continuous if its corresponding ISPES is 
continuous. 

If a SPES is not continuous, one can identify a couple of neighboring states that surround 
a discontinuity of the mean-value of one of the possible operators on the corresponding ISPES, 
thus localizing this discontinuity on the SPES. 

For a given discontinuity, at least one of the mean-values of the usual multipole moments 
must be discontinuous. We call the lowest-order discontinuous multipole moment operator the 
signature of the discontinuity. 

2.2. Density distance 

As a way to numerically estimate the difl'erence between two neighboring states and 
we define the density distance Dpp' as 



Ve,^,Vxe7,V;e{l,A^,}, 
'^M^iqxj + d)\QAiiW{qx,j + d)) 

= <i/''(?x,>)|(2^^|iA'(?x,j)) 



(8) 



For the ISPES of the previous example, the continuity condition is 




= {4''iq2o(x))\QAn\if''(q2o(x))} 



(9) 




(10) 
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with p(r) and p'{r) being the local spatial density of the |t/r) and |(/r') states, respectively. This 
distance is dimensionless, and can be seen as the geometrical difference between two states, 
measured in nucleons. Any state of A nucleons has a density distance of A with the void state 
( / dT^p(r) = A and p'(r) = 0). 

2.3. Maximum density distance 

We define the maximum density distance ZJtnax of a SPES as 

£>max = max{Dp.p) 

with i e {I, N],j e Ni, (11) 

A^, being the set of closest neighbors of states li/r),-. In the case of a regular mesh, two points 
are closest neighbors if only one of their coordinates differ and if this difference is minimal. In 
the case of a non-regular mesh, two points are considered as closest neighbors if there exists a 
Delaunay cell containing them |24]. Being given a A^^. -dimensional set of points and an Euclidean 
distance function, the Delaunay triangulation is unique (if no points are co-circular). The distance 
function used for the Delaunay cells computation has to be carefully chosen in order to take into 
account variations in all Nc possible directions. For example, to find the closest neighbors in the 
case of a non-regular 4-dimensional mesh with constraints on (Q20,Q30,Q40 and Qeo), we have 
used the following pseudo-Euclidean distance function: 

d' = + 



lb / \ lb^/2 

I IbM I Ib^ I ^ ^ 

with the values of q2o, Qio, q40 and ^60 being given in b, b-'''^, b^ and b^, respectively. A 2- 
dimensional example of a Delaunay triangulation using the usual Euchdean distance is given on 
Fig. [3] 

An example of a discontinuity for the symmetric fission of ^^^Fm is shown on Fig. |4] The 
evolution of the HFB energy is rather smooth (top panel of Fig|4|i. There is no way to suspect a 
possible discontinuity just by considering the HFB energy curve. However, a spike of the Dpp' 
value (bottom panel of Fig. |4]l is observed for q2o - 260 b and ^20 = 261 b, consequence of a 
discontinuity. We can see that the mean value of the hexadecapole operator Q40 changes at this 
discontinuity (middle panel of Fig. |4]l, so we can conclude that this discontinuity has Q^q for 
signature. Since the evolution of the energy is monotone, no saddle point is visible, but there 
has to be a saddle point between both valleys in the Q40 subspace. This is an example of a 
missing saddle point. The local densities before and after the discontinuity are shown on Fig. |5] 
The main visible change in the local densities seems to be a stronger neck between the nascent 
fission fragments. 

An other example of a discontinuity is shown on Fig. |6] for the symmetric first barrier of 
^'^''Pu. We would like to stress that this axial 1 -dimensional SPES is shown as an example only, 
and must not be considered as a prediction of the real first barrier of ^"^^Pu. The influence of the 
triaxial degree of freedom is neglected here, and it has been shown that this degree of freedom 
significantly lowers the first barrier of several actinides Jls, 26, 27, 28, 29, 3^ 31, 32]. The 



situation is very similar to the one on Fig. |4| the evolution of the HFB energy is smooth (top 
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Figure 3: Example of a Delaunay triangulation in a 2-dimensional space, using the usual Euclidean distance. Points E, 
F and L belong to the same Delaunay cell. In a W-dimensional space, a set of (W + 1 ) points is a Delaunay cell if its 
circumscribed circle/sphere/. . . contains no other point. The circumscribed circle of the Delaunay cell EFL is shown in 
red. 

panel of Fig. |6|, but a spike of the Dpp' value (bottom panel of Fig. |6]l indicates that there is 
a discontinuity between q2o = 48 b and q2o = 50 b. The signature of this discontinuity is the 
hexadecapole operator Q^q, since its mean value is discontinuous as shown in the middle panel 
of Fig. |6] The local densities corresponding to the states before and after the discontinuity are 
shown on Fig.|7] The visible difference between these states seems to be the quantity of matter 
around the z - plane. Since the discontinuity occurs at a saddle point, this point is what we 
have called a fake saddle point in the Introduction. 

Since the step for the ^20 constraint value is smaller for the SPES of Fig.|4]compared to the 
SPES of Fig. |6l it is important to remark that the average value of Dpp> for continuous points is 
smaller (around 0.3 for ^^''Fm versus 1.0 for ^"^''Pu). To check the continuity of a SPES, the used 
steps must be small enough to be able to see spikes corresponding to discontinuities. 

3. Obtaining a continuous SPES 

As far as we know, all of the discontinuities we have had to deal with led at least to a density 
distance of several nucleons. It is then safe to consider that a SPES is continuous if we have 
Dmax < 2. In order to reach this level, all of the discontinuities have to be "cleaned". There are 
two ways to do so. 

3.1. Additional dimension 

The first way to remove a discontinuity is simply to globally extend the mesh with the sig- 
nature of this discontinuity. For example, if a given one-dimensional SPES along Q20 has a 
discontinuity with ^40 for signature, the two-dimensional SPES in the {Q20, Q40) subspace will 
not have any discontinuity with Q40 for signature. This method is brute-force, and can lead to a 
high number of points. As an example, we have extended the SPES shown on Fig.|6]by taking 
the signature of the discontinuity Q40 as a constraint. The resulting SPES is shown on Fig. [8] On 
this figure, the topology of the two valleys (corresponding to the first and second well) is very 
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<Q,„> [b] 



Figure 4: Example of a discontinuity for the symmetric scission of^sfipm. The upper panel shows the total binding 
energy, the center panel shows the mean value of the hexadecapole operator Q40, and the lower panel shows the maximum 
total density distance between each point and its closest neighbors. The cylindrical symmetry is conserved. 



similar to the one of the schematic valleys shown on Fig.|2] One can see that the automatic basis 
optimization can be a source of rather small local maxima of Dpp> values, since these maxima 
values are lowered when disabling the basis optimization. Anyway, these local maxima are small 
enough to let us consider that the SPES shown on Fig.[8]are continuous. 

3.2. Connecting points 

The second way to remove a discontinuity is to locally extend the mesh with some missing 
points at a discontinuity. For the schematic discontinuity of Fig.|2l one could simply add to the 
SPES the connecting points between valleys A and B of the panel (f). To find these connecting 
points, one calculates a small two-dimensional SPES large enough to include points a and b 
like the one plotted in panel (c), one finds the saddle point between the valleys by calculating 
the discrete local derivatives for each point, or by a more subtle technique if needed, like for 



example the immersion technique BI8L 13311 . and one calculates the least-energy paths from this 
saddle point to the bottom of each of the valleys. The immersion technique is a fast and reliable 
method to find the saddle point between two given points on a regular mesh. Its principle is 
to mimic the way a liquid would fill a valley before flowing into an other one through a saddle 
point. 

The schematic result is shown on Fig.|9] It is worth mentioning that (22o) is no longer a good 
variable for a one-dimensional description, for two reasons. Firstly, it is possible that the least- 
energy path "backbends", i.e. there can be multiple points with the same {Qio) value. Secondly, 
the connecting points are not anymore minimizing the total binding energy for a given {Qio) 
value. This is why a curvilinear abscissa is used instead of {Qio), and why the original points 
between b and c are discarded and replaced with points from the least action paths. The resulting 
set of points can no longer be called a SPES, but can be called an effective PES. 
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z [fm] 



Figure 5: Local densities for -^*Fm in the (z, r) plane before (top) and after (bottom) the discontinuity shown on Fig.|4] 
The isoHnes are separated by 0.02 fm"'. 

This method is rather convenient in the case of one-dimensional SPES, but can be hard to 
automatize for SPES of higher dimension. The fact that the constraint variables are no longer 
good variables can also be problematic when trying to perform dynamical propagations on the 
resulting effective PES. 

5.5. Discussion 

We would like to stress that for a given system, the dimension of the smallest continuous 
SPES is the dimension of the smallest corresponding continuous ISPES. If the continuous de- 
scription of a physical process needs A^^ degrees of freedom, the dimension of a corresponding 
continuous SPES must obey N > Np. Trying to obtain a continuous description with an insuf- 
ficient number of degrees of freedom with a change of variables like in the "connecting points" 
method presented above can be dangerous for the physical interpretation of the resulting quan- 
tities like the inertia tensor, mainly because the fact that some points do not minimize the total 
binding energy anymore for a given (22()) value forbids the use of this latter as a good variable. 

Using the "additional dimension" method presented above does not lead to such problems, 
but can be numerically difficult. For some physical processes, having = Np can mean that 
the number of points of the continuous SPES is huge if the considered mesh is rectangular and 
regular As a way to lower this number of points, one can restrict the SPES to the solutions with a 
low energy (non-rectangular mesh), and/or use multidimensional adaptive meshing (non-regular 
mesh). 

Conclusion 

Discontinuities are not easily visible, and have important consequences. Any result extracted 
from a discontinuous SPES should be considered with extreme caution. The method to obtain 
Dma\ that we have presented in this paper is simple, and allows to check with good confidence 
if a given SPES is continuous or not, without any additional self-consistent calculation. The 
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Figure 6: Example of a discontinuity for the symmetric first barrier of ~'"'Pu. The upper panel shows the total binding 
energy, the center panel shows the mean value of the hexadecapole operator Q40, and the lower panel shows the maximum 
total density distance between each point and its closest neighbors. The cylindrical symmetry is conserved. 



value of Dmax should be specified for any SPES, like the conserved symmetries and the level of 
convergence reached. 

In the near future, we plan to perform dynamical calculations on SPES for the description 
of the fission process and the prediction of some fission fragment properties, using a Time- 
Dependent Generator Coordinate Method |34, 35tl, and the continuity of the produced SPES 



will be carefully checked using the method presented here. One of the main advantages of this 
method is that it can be fully automatized, allowing an on-the-fly detection of discontinuities 
during the production of a SPES. 

In this study, the definition of Dpp> has been chosen for practical reasons (mainly its numerical 
simplicity), but other expressions can be used to estimate the "distance" between two neighboring 
states, involving for example the non-local density 

0;^, = J Jd}'di'Wy)-p'(r'y)\ (13) 

or the overlap between neighboring states 

Such expressions can reveal a discontinuity in non-local fields, but their computational cost may 
be important. One should also note that the proposed expression does not take into account time- 
odd quantities. When calculating a SPES for an odd nucleus, another expression should be used if 
one is looking for discontinuities in time-odd quantities. In Ref. lUsll . the authors state that "[. . . ] 
minimization does not necessarily lead to the wrong solution; in many simple problems it will 
give the correct saddle point. Unfortunately, one cannot know a priori whether a true or incorrect 
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Figure 7: Local densities for ^''"Pu in the (z, r) plane before (top) and after (bottom) the discontinuity shown on Fig.|6] 
The isolines are separated by 0.02 fm"'. 



solution will be found". Thanks to the simple method presented in this paper to calculate Dmax, 
it is now easy to estimate if and where a discontinuity has occurred while calculating a SPES. 
One can consider that a SPES with a low enough Z)max value does not contain such "incorrect 
solutions" (fake or missing saddle points), and that the "Problem SPES" from Introduction is 
fixed. By comparison, there is to our knowledge no numerical method to fix "Problem NSPES" 
from the introduction, in other words, to estimate if partial local extrema obtained by non-self- 
consistent methods in a given subspace remain partial local extrema in bigger subspaces. 
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figures. 
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Figure 9: Illustration of the "connecting points" method. The saddle point a is found, then points alo b and a to c are 
found as least action paths. Original points before b, original points after c, and both least action paths (a to b and a to 
c) are then joined to form a continuous path. The constraint value ^20 from Fig.|2]has been replaced with a curvilinear 
abscissa i in lower panel. 
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